Spin transverse force and intrinsic quantum transverse transport 
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The spin-orbit coupling may generate spin transverse force on moving electron spin, which gives 
a heuristic picture for the quantum transverse transport of electron. A relation between the spin 
and anomalous Hall conductance and spin force was established, and applied to several systems. It 
was predicted that the sign change of anomalous Hall conductance can occur in diluted magnetic 
semiconductors of narrow band and can be applied to identify intrinsic mechanism experimentally. 
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PACS numbers: 72.25.Dc, 72.20.-i, 72.10.-d, 85.75.-d 



I. INTRODUCTION 

The theory of anomalous Hall effect has a long history 
since 1950si It was realized that several different mech- 
anisms are contributed to the total Hall conductance in 
ferromagnetic metals and semiconductors. Except for 
the skew scattering and the side jump from the impu- 
rity scatteringjSiiii^ an intrinsic mechanism tells that the 
spin-orbit coupling in the electronic band structure of 
the system may induce a non-zero Berry phase or mag- 
netic monopole in the momentum space and cannot be 
neglected in the transverse transport of electrons, espe- 
cially in the diluted magnetic semiconductors*^ On the 
other hand, the spin aspect of transverse transport was 
also studied recently, in which an external electric field 
may drive electrons to form a transverse spin current in 
the systems with spin-orbit coupling even in paramag- 
netic electronic systems .SaLiS These two effects reflect 
the charge and spin aspects of electron transport, respec- 
tively, and have some common features as their physical 
origin stems from the same spin-orbit coupling of con- 
duction electrons. More and more experiments support 
the intrinsic mechanism^iiiSii^ Since the spin-orbit cou- 
pling reflects the interaction between the electron spin, 
momentum and external or environmental potential, it 
was shown that this coupling, as an extension of Ehren- 
fest's theorem in quantum mechanics, may generate a 
spin transverse force on the spin current instead of the 
Lorentz force on electric current in a magnetic field in 
conventional Hall effect ^ This spin transverse force pro- 
vides a new route to get insight of the transverse motion 
of electrons. 

In this paper we focus on the relation between the spin 
force and intrinsic quantum transverse transport. For the 
electronic system with spin-orbit coupling and exchange 
coupling, it was found that a spin transverse force exerts 
on moving electron spins in the anomalous transverse ef- 
fect, as a result, a set of formula based on spin force were 
derived for intrinsic spin and anomalous Hall conduc- 



tance. As an intrinsic feature it was predicted that the 
sign of the anomalous Hall conductance can be changed 
when the system breaks both the bulk and structural in- 
version symmetry and time reversal symmetry. It is also 
shown that spin Hall effect in the two-dimensional (2D) 
Luttinger system is equivalent to the intrinsic anomalous 
and spin Hall problems in a ferromagnetic metal, and 
its robustness against impurities can be also understood 
very well from the point of view of spin force. 



II. GENERAL FORMULA 

To develop a general formula for spin and anomalous 
Hall conductance, we start with an effective Hamiltonian 
for electrons with spin 1/2, 



H = e(p) + ^2 d a (p)a a , 



(1) 



where e(p) = p 2 /2m* is the kinetic energy with the 
band electron effective mass m* and a a are the Pauli 
matrices. d a (p) are the the momentum-dependent co- 
efficients which describes the spin-orbit interactions and 
exchange interaction of magnetic impurities. The energy 
eigenvalues are E pn = e(p) + pd(p) with ji = ±, and 

d = y/ a* + d 2 + d 2 , and the corresponding states are 

denoted by \p,/J,). Using the Heisenberg equation, the 
kinetic velocity operator (j = x,y,z) is defined as 



Pi 



ih m 



E 

a=x,y,z 



dd a 
dPj' 



(2) 



where the first term is the canonical velocity, and second 
part can be regarded as the spin gauge field or anomalous 
velocity, (e/m*c)Aj = Y, a =x,y,z®da/dPj<r a - As an ex- 
tension of the Ehrenfest's theorem^ the spin-dependent 
force is introduced as the derivative of kinetic momentum 
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operator with the respect to time, 

m* 2m* dd a 

F j = —[v j ,H] = —e aPl —d p a 1 . (3) 

This spin force is a purely quantum quantity, and has no 
classical counterpart. 

The quantum transverse transport of electrons caused 
by a weak electric field E can be calculated within the 
framework of linear response theory. The Kubo formula 
for the d.c. conductivity gives^ 



E 



Uw - fw') Im ((pmI v i W) (pm'I v j W)) 



(E pfl - Eppi) (E, 



E 



vp 



iS+) 



(4) 

with 5 + — ► + , and the Dirac- Fermi distribution function 
f p4l = 1/ {exp[/3 (E^p) - n)] + 1} . With the formula of 
the kinetic velocity in Eq.@ 



a=x,y,z 3 



(5) 



for /.t 7^ // and furthermore an identity can be proved for 
this system, 



Im((p/x| cr Q \pp!) (pfj,'\ ap \pfi)) = Mea/3 7 -7-> 



(6) 



for /i 7^ //. We limit our discussion to the case that 
two spectra are non-degenerated in the whole momen- 
tum space, i.e. , d in the whole momentum p space 
such that 5 + can be taken to be zero before the integral 
over pi& Thus the intrinsic transverse conductance can 
be expressed as 



e 2 fi (fp-- ^ fp-+) ddadd^ 



2tt ^ 

p 



d 3 " api d Pl d Pj ^- 



(7) 



This recovers the conductance formula in terms of the 
Berry curvature ^Li&iS On the other hand, the anticom- 
mutators of the spin force Fj and spin gauge field Ai 
gives 



Tr[{F i! A}] = 



tic 



dd a dd 



dpi dpj 



(8) 



where the trace runs through the spin variables, 
{PlA ' ' ' IPA 4 )- I n this way we established a relation 
between the electric conductance and spin force, 



Onsager relation. It is noticed that Tr[{Fj, Ai}] = if 
any one of d a is zero. 
The Kubo formula for spin Hall conductance is written 

as 

7 eft (Jpp - f w >) Im ( Pf i\ J'l \pp!) (p/i'l Vj \pfi) 



E„ 



( 10 ) 

where spin current operator J 4 7 is defined conventionally 
as J? = {h/A) {vi, (T 7 }, note that this choice a natural 
one but not a unique one in the presence of spin-orbit 
coupling since these is no continuity equation for spin 
density as is the case of for charge density. Following 
the calculation mentioned above, from Ea. (|lU|l one can 
obtain^ 



t _ eft2 \ ' Up- ~ fp,+) de ddp 



(ii) 



Similarly, the spin conductance can be also expressed in 
terms of the spin force, 



lQm*Q 



E 



Up - — fp,+) 

d 3 



Ti[{Fj,J?}], (12) 



which is given by the anticommutator of the spin force 
and spin current operators J? with 



(13) 



On the other hand, due to the spin-orbit coupling, an 
electric field can also induce a non-zero spin polarization. 
Its linear response to the field is, 



2 - K)ij=o) l E 3 



, h 3 



32m*tt ^ 

p 



Up - — fp,+) 
d 3 



Tr[{F j)( 7 7 }].(14) 



In this way we have established an explicit relation be- 
tween the spin force and the quantum spin and charge 
transverse transport in a system with spin-orbit coupling. 

In general, for any observable O, its linear response to 
an external electric field E, (Oi) = \ijEj, where, from 
the Kubo formula, 



Xij 



16m*n 



16to* 2 c£1 



E 



Up,- — fp,+) 

d 3 



Tr [{Fj,Ai}}. (9) 



III. APPLICATIONS OF THE GENERAL 
FORMULA 



Since Tr[{Fj,Ui}] = ^Tr[{Fj, Ai}] , the anomalous 
Hall conductance is determined by both the spin force 
and the spin gauge field, or anomalous part of the veloc- 
ity. This fact reflects the physical origin of spin-orbit 
coupling in this effect. Note that Tr[{i^-,^}] = — 
Trfl-Fj, Aj}], so <xy = — <Xj, for i ^ j, that is the so-called 



A. 2D ferromagnetic system with the Wurtzite and 
zinc-blende structures 



Now we apply the general formula to several systems 
in semiconductors. We first consider an effective Hamil- 
tonian for a two-dimensional ferromagnetic system with 
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the Wurtzite and zinc-blende structures, 



with 



H 



2m 



- + aof(p){p y cr x - p x (T y ) + h a z 



(16) 



The exchange field due to the magnetic impurities or 
Coulomb interaction is taken to be uniform, as an mean 
field effect, and is normal to the plain. f(p) = 1 for the 
Wurtzite structure, and f(p) = p 2 for the zinc-blende 
structure such as Hgi-^Mn^Tei^S The components of 
spin gauge field A in this system are 



A> 



Ay 



m cclq 



m cclq 



df (., df 



(17) 



(18) 



and the spin force is, 



4to* 



alfJl x z 



2m*aoho 



V p (fp-a). (19) 



From the formula we notice that a spin transverse force 
arises due to the spin-orbit coupling if there exists a spin 
current along the electric field and also an in-plain spin 
force arises due to the interference between the spin-orbit 
coupling and the exchange field. To calculate the intrinsic 
Hall conductance, we notice that the anticommutator, 



Tr [{F y ,A x }} 



8m* 2 ch 
he 



2 of[d(pf)/d P }, (20) 



which vanishes when ho = 0. Thus if the system does not 
violate the time reversal symmetry no anomalous Hall 
current circulates. Assuming at zero temperature and 
taking the bottom subband to be occupied and the top 
one to be empty. The anomalous Hall conductance is 



p F = V2{[m* 2 (4m* 2 + hi + 2alfim*)] 
+m* (a 2 m* +^)} 1/2 . 

For the zinc-blende structure, i.e., f(p) = p 2 , 

z _ eh 2 (h -m) 2 



1/2 



where 



xy 32to*7t ag?72 

1 1.1 i 4 2\ 1 / 2 

m = (h +P F a ) 



(25) 



(26) 



(27) 



with 



p F = {2m*( M + [hi - Aalm* 2 (h 2 - M 2 )] 1/2 ) 

/(l-4aX 2 )} 1/2 - (28) 



B. 2D ferromagnetic electrons gas with both 
Rashba and Dresselhaus coupling 

In some materials the spin orbit coupling has two dif- 
ferent contributions when a system has the structural 
and bulk inversion asymmetry. Here we consider a 
two-dimensional ferromagnetic electrons gas with both 
Rashba and Dresselhaus coupling, 

p 2 

H = - \{p x <Jy -PyO x ) ~ P(p x a x -PyCJy) + h (J z . 

(29) 



In the language of d a in Eq.QJ, 

d x = Xp y - (3p Xi 



(30) 



e 

2h 



(l - h / [a 2 x 2 + hi 



2lV2 



(21) 



where x 2 — p 2 F j 2 at the Fermi surface. If ooPf/ >> 
ho, the anomalous Hall conductance is almost quantized, 

2 

a xy ~ — m-. This is consistent with the theory of Berry 
curvature m& 

According to the Ea. (|12fl . one can obtain the intrinsic 
spin Hall conductance 



eh 2 x - 

4m*tt 

p 



Up,- - f P ,+)P 2 x a 2 of 2 



[p 2 alf 2 



hlf 2 



(22) 



For the Wurtzite structure, i.e., f(p) = 1, assuming the 
Fermi energy /x lies in the gap, the intrinsic spin Hall 
conductance is 



where 



z eh 2 {hp-rjiY 
J ' xy ~ 16to*tt alni 

Vi = (h Q + PF a o) 



(23) 



(24) 



dy = ~\p X + 



d z = h . 



It is straightforward that 



Tr[{F v ,A x }} 



8m* ch 
he 



(A 2 -/? 2 ), 



(31) 



(32) 



(33) 



which is independent of p. Using the formula for the Hall 
conductance, we conclude immediately that the anoma- 
lous Hall conductance vanishes at ho = 0. The Hall con- 
ductance will change its sign near the point of X 2 — (3 2 . 
Since /io / 0, the band spectra opens a finite gap at 
p = 0, and degeneracy of two electron bands is removed. 
If one assume that the Fermi surface is below the the up- 
band of the spectrum, the anomalous Hall conductance 
<j xy reads 



2h 



n (A 2 - d 2 ) (1 - 7 ) , 



(34) 
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where 

1 f 2 * |A 2 -/? 2 | , , 

1=7^ deb ! — , 35 

2W0 n(<p)[i+ p* F (<b)n(<t>)/h 2 ] y2 

with 

II ((b) = X 2 + f3 2 — 2X(3 sin 20, (36) 

and pf is the Fermi momentum, 

p F = V2m*{p + m*XP (A/3 - 2 sin 2(b) + [hi - 

+ (m*\ 2 (3 2 +n- 2m*A/3sin20) 2 ] 1/2 } 1/2 . (37) 

The dependence a xy on (3 of numerical calculations is 
plotted in Fig.l. The parameters are taken as 2 ^: A = 23 
meV nm/fi, ho = 1-38 meV, /i = 1.38 meV, m* = 0.9too, 
and /3 is comparable with A. It is shown that the anoma- 
lous Hall conductance will change its sign near the point 
of A 2 = /3 2 . If pf |A — /3| » /loi o"xy is approximately 
equal to — sgn(A 2 — /3 2 ) 

If one assume that the Fermi surface is below the the 
up-band of the spectrum, the intrinsic spin Hall conduc- 
tance can be also obtained 

a* =^J— f 2 * {\ 2 -f3 2 )(h a - V ((h)) 2 cos 2 (j> 
° XV 4m(27r) 2 y o II ((b) 2 7] {(b) 

.(38) 

While the Fermi energy is above the gap, the intrinsic 
spin Hall conductance is 

z _ fi 2 e^_ f 2v (A 2 - P 2 ) cos 2 (b 
° xy ~ 4m(2n) 2 Jo n(0)V?/ + 

x (h 2 - v _ r,+)(r, + -r?_), (39) 

with 

ij± = [h 2 Q + p 2 ± n((b)} 1/2 , (40) 

and 

p± = V2m*{n + m* A/3 (A/3 - 2 sin 20) =F [hi - M 2 

+ (m*A 2 /? 2 + /x - 2m*A/3sin20) 2 ] 1 / 2 } 1 / 2 . (41) 

The numerical calculation of the intrinsic spin Hall con- 
ductance a z xy with (3 is plotted by the dashed line in Fig.l. 
When sgn(A 2 — f3 2 ) = +1, the intrinsic spin Hall conduc- 
tance a% y is positive; sgn(A 2 — /3 2 ) = —1, a* y is negative. 

The sign change of anomalous Hall conductance as 
that of spin Hall conductance^ may be observed experi- 
mentally in some diluted magnetic GaAs quantum wells 
where the Rashba and Dresselhaus coupling are usually 
of the same order of magnitude. The Rashba coupling is 
adjustable by a gate field perpendicular to the electron 
gas plane. Thus there is no practical difficulty to achieve 
the situation near A = ±/3. Since the Rashba coupling 
has a relation to the gate voltage, A = Aq + SXV g , the 




P 

FIG. 1: Variations of the conductivity a xy and the spin con- 
ductivity a% y with the coupling parameter (3 (meV nm/h). 
The parameters are A = 23 meV nm/fi, ho = 1.38 meV, 
fj, — 1.38 meV, and m* = 0.9mo. The solid line corresponds 
to (T X y, and the dashed line to a% y . 

gate field can be used to control the direction of intrin- 
sic anomalous Hall current. If the extrinsic contribution 
to the Hall conductance is comparable with the intrinsic 
one, at least a jump of the anomalous Hall conductance 
should be observed even if the Hall current could not 
change its direction. This effect can be identified as the 
intrinsic mechanism for anomalous Hall effect. 

C. 2D Luttinger model 

We now turn our attention to the Luttinger Hamil- 
tonian for spin 5 = 3/2 holes in the valence band of 
centrosymmetric cubic semiconductors^ 

H=^-( ll + ^ 2 )p 2 -^(p-S) 2 , (42) 
2m \ 2 J m 

where 71, 72 are material-dependent parameters and m 
is the electron mass. In terms of SO (5) Clifford algebra, 
the Hamiltonian (|42|l can be cast into^ 

H = e(p) + d a T a , (43) 

where e (p) = 7ip 2 /2m, r a (a = 1, 2, • • • , 5) the five 
Dirac T matrices and d a the five d-wave combinations 
of p are given in Ref^. In the two-dimensional case, 
for the first heavy- and light-hole bands, the confinement 
in a well of thickness a is approximated by the relation 
(p z ) = 0, (jp 2 ) ~ (irh/a) 2 ~ In this case, c?i = d 2 = 
and r a (a = 3, 4, 5) forms reducible representation of an 
SO (3) Clifford sub-algebra»i£. In the following, one 
introduces a new representation under which the ex- 
pressions of S matrices and r a (a = 3, 4, 5) are given 
in Appendix. Under this new representation, the two- 
dimension Luttinger model is block diagonal 
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where 



= e (p) + /u d a a a , 



a—x,y,z 



with /i = ±1, and 



e 0) = 7T- fal + Py) > 
2m yj 



, V372 

^a; — PxPy-i 



, ^72 / 2 2\ 



(45) 

(46) 
(47) 
(48) 



o 






72 
m 



(49) 



(see Refsi 2 ^ 2 ^* 2 ^) These two effective Hamiltonians are 
connected by the time reversal operator 6 = ~ia y K (K 
is the complex conjugate operator that forms the com- 
plex conjugate of any coefficient that multiplies a ket): 
H+ = G.ff_0 _1 . Thus the energy eigenstates of H are 
at least doubly degenerated. has the same form as 
that for the ferromagnetic metal with exchange coupling 
since d z (p) in contains a non-zero term, — ^ (Pz)> 
but the Pauli operators <r a here do not represent a real 
spin. To calculate the linear response of the spin current 
J? = I {dH/dpi, S z } to an electric current Vj = dH/dpj, 
the Kubo formula gives 



■xy 



M=±l 



Au xy ' 2e 



(50) 



where crl^ and a 



'xy 



conductance of H^, respectively. 



are the anomalous and spin Hall 
The anomalous Hall 



conductance cr.y = X^=±i °xy - These relations are valid 
even if the non-magnetic disorder and interaction are 
taken into account. It was known that the spin Hall con- 



ductance is invariant o~xy +1 ^ 



but the anomalous 



Hall conductance changes its sign a. 



(+1) 

xy 



-a. 



(-1) 

xy 



under 



time reversal. Therefore the calculation of spin Hall con- 
ductance in the 2D Luttinger model is reduced to those of 
spin Hall conductance and anomalous Hall conductance 
in a ferromagnetic metal or band insulator of H^. For a 
numerical calculation the material-specific parameters of 
band structure for GaAs are adopted as 71 = 6.92 and 
72 = 2.1 the thickness of quantum well a = 8.3 nm cor- 
responds to the gap between the light- and heavy-hole 
bands at the T-point AE = 40 meV. Using the formula 
of Eq. @ and © &xy 1 @xy 

, and a z xy are calcu- 
lated numerically in terms of chemical potential /1 as 
shown in Fig. 2. It is noted that, for an infinite confine- 
ment, {pi) —> +00, the eigenstates of if M becomes fully 
saturatediSi As a result the spin-orbit coupling is sup- 
pressed completely, and both spin and anomalous Hall 
effect vanish, which is in agreement with Bernevig and 
Zhang^i 



FIG. 2: Variations of ai^, a. 



and with Fermi level 



pu (meV). The parameters are a = 8.3 nm, 71 = 6.92 and 
72 = 2.1. The solid line corresponds to <r% y , the dotted line 

to <Txi +1 \ and the dashed line to crit 1 '. 



IV. DISCUSSION 

Here we may present a heuristic picture for electronic 
transverse transports in these systems and the effect of 
disorder from the point of view of spin force. The effect 
of disorder on anomalous Hall conductance has been in- 
vestigated in RefsiS&SS. For simplicity we focus on the 
system with the Wurtzite structure, i.e., / = 1 in Ea. l|16() 
or j3 = in Ea. H29(l . The spin-orbit coupling induces the 
spin force, which contains two parts: the transverse force 
on moving electron spin, 



4m* 2 A 2 



ft 2 



-Jf x z, 



(51) 



and the exchange coupling interacting also induces a 
spin force within the plane, which is relevant to the spin 
polarization, 



'2 = 



a x x + a v y] . 



(52) 



If the disorder potential Vdisorder is taken into account, 
in a steady state, the spin force may reach at balance, 



1 

ih 



-A, H + V d 



isorder 



= (Fi + F 2 )=0. 



(53) 



This result is independent of the non-magnetic disorder 
and interaction because the spin gauge field commutes 
with non-magnetic potential Vdisorder- From the spin 
force balance we have a relation between spin current 
and spin polarization, 



(Jx) 

(Jy) 



hhQ 

2m* A 
hho 
2m* A 



(54) 
(55) 



It becomes obvious that the spin Hall current vanishes 
in the case of ho = 0, which is consistent with previous 
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complicated calculationsi^SiSiiSiS Of course the purely 
intrinsic responses of spin current and spin polarization 
without impurities do not satisfy this relation, and the 
extrinsic contributions have to be included to reach the 
balance. Thus, this gives a clear picture for anomalous 
Hall effect in ferromagnetic metal: When the external 
electric field is applied along the x axis, it will circulate 
an electric current J c , x , and also a spin current J x since 
the charge carriers are partially polarized. The spin-orbit 
coupling exerts a spin transverse force on the spin cur- 
rent, JJ, and generate a drift velocity or the anomalous 
Hall current J c ,y From the Rashba coupling the spin 
polarization tends to be normal to the momentum or 
electric current. The electric current J c ^ x along E in- 
duces a non-zero (a y ) and the anomalous Hall current 
J CiV induces non-zero (<j x ) . These non-zero spin polar- 
ization maintains the balance of spin transverse force, 
and further a non-zero spin current in a steady state. 
Thus the anomalous electronic transverse transport is ro- 
bust against the disorder in the ferromagnetic metals and 
semiconductors. This picture can also be applied to un- 
derstand the spin Hall effect in the quasi-2D Luttinger 
system. 
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Now the Dirac matrices T a (a = 3, 4, 5) become 



r 3 = — (S X Sy + SyS X ) =<J z (g)(T x 



(A3) 



r 4 = — (S z x - S z y ) = <J z (g a 



(A4) 



(A5) 
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Appendix 



In order to write the two-dimensional Hamiltonian in 
the form of Eq. (|44|l . one introduces a new representation 
under which S matrices are expressed as 



r 34 = 7r7[r 3 ,r 4 ] = i®<7„ 



(A8) 



S z 



f-l o o 0\ 

\ 
-\ 

V o o o 1 1 



(Al) 



s z — — r 3 4 — 2^12- 



(A9) 
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